Thermodynamics of rotating Bose gases in a trap 
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Novel ground state properties of rotating Bose gases have been intensively studied in the context 
of neutral cold atoms. We investigate the rotating Bose gas in a trap from a thermodynamic 
perspective, taking the charged ideal Bose gas in magnetic field (which is equivalent to a neutral gas 
in a synthetic magnetic field) as an example. It is indicated that the Bose-Einstein condensation 
temperature is irrelevant to the magnetic field, conflicting with established intuition that the critical 
temperature decreases with the field increasing. The specific heat and Landau diamagnetization also 
exhibit intriguing behaviors. In contrast, we demonstrate that the condensation temperature for 
neutral Bose gases in a rotating frame drops to zero in the fast rotation limit, signaling a non- 
condensed quantum phase in the ground state. 

PACS numbers: 03.75.Hh, 03.75.Lm, 05.30.Jp, 75.20.-g 



The investigation of rotating quantum gases is one of 
the central topics in the study of superfluidity and su- 
perconductivity. So far, research attention in this re- 
gard mainly focuses on the physics of quantized vortices 
in superfluids or superconducting samples. Especially 
in recent years, several experiments using rotating Bose- 
Einstcin condensates (BEcs) of trapped alkali atoms have 
provided spectacular illustrations of the notion of quan- 
tized vortices [lH6| , and thus stimulate enormous interest 
in the properties of rotating condensates Q • In these ex- 
periments, the atomic gas is equivalent to being confined 
in a rotating frame, where the Hamiltonian of a single 
particle is given by 



quantum Hall regime) [9-41 1| to the highly correlated non- 
condensate state (referred as the fractional quantum Hall 
regime) [IH, [l3[ takes place. 

An alternative approach to spin up neutral atomic 
gases is introducing a "synthetic" magnetic field (l^4l6| , 
instead of rotating the frame. This approach creates ef- 
fective gauge potentials A for atoms by dressing them in 
a spatially dependent manner with optical field that cou- 
ples different atomic internal states. The effective Hamil- 
tonian for a trapped neutral atom in the "synthetic" field 
is like that of charged particles in the magnetic field, 
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with wo denoting the frequency of harmonic potential in 
the x,y plane and r 2 = x 2 +y 2 . ft is the rotational angular 
frequency around the z axis, which results in an effective 
vector gauge potential, A = (—MCly,MD,x,0). As the 
rotation frequency fl increases, one or several vortices are 
formed in the condensate. With more rapid rotation, the 
Abrikosov vortex lattice can be observed 0, S @ • 

According to the Hamiltonian (JXJ) , the regime where 
fi near the trap frequency ujq (called the fast-rotation 
regime) is of special interest 7, 8). As ft -> u)q, the single- 
particle energy levels become macroscopically degenerate 
and the rotating dilute gas is expected to exhibit novel 
quantum phases analogous to the quantum Hall state of 
electrons in the strong perpendicular magnetic field, if 
interactions between atoms are taken into account |9h13|] . 
During this process, a quantum phase transition from the 
supcrfluid state (Bose-condensed state at the mean-field 
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Vf (P - A) + -MwZ(x* + y 1 ) + Kxtra(r). (2) 

Here T4 x t ra is the modification to the trapping poten- 
tial, which is related to the optical field where A is 
induced [l6j]. Both A and V cxt ia can be engineered 
flexibly through constructing the position-dependence of 
the dressed atomic states [la, LL2I- Quite recently, the 
light-induced magnetic field is realized in experiments, 
with producing either uniform vector potential 17| or 
spatially- varied one [lSj. The latter corresponds to a 
nonzero magnetic field which can generate vortices in 
the ground state of BEcs. and this phenomenon has just 
been indeed observed [18| . An advantage of this optical 
approach over rotating gases is that the synthetic field 
can be significantly large, making possible the study of 
the fractional quantum Hall regime. 

Both approaches mentioned above acquires a vector 
potential for neutral atoms, so as to imitate charged par- 
ticles in the magnetic field. Actually, the charged Bose 
gas (CBG) has already been studied for more than half a 
century. Earlier in 1950s, it showed that the ideal CBG 
exhibits essential equilibrium features of a superconduc- 
tor |l9j, l20|. Furthermore, it was indicated that an arbi- 
trarily small value of the magnetic field can eliminate 
Bose-Einstein condensation (BEC) in a 3-dimensional 
ideal CBG, while the orbital motion results in extremely 
large Landau diamagnetism and even leads to Meissner 
effect at low temperatures [2£}-[2jj. 
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In this paper, we study properties of rotating Bose 
gases in a trap from a thermodynamic perspective. With- 
out loss of generality, we consider a ideal CBG trapped 
in the harmonic potential as an example. The obtained 
results are applicable for neutral Bose gases either in 
a rotating frame or in a synthetic magnetic field. We 
demonstrate that these two cases display distinct ther- 
modynamic properties. 

The ideal Bose gas model. - Let us consider an ideal 
Bose gas with mass M and charge q. The whole system 
is placed in a uniform magnetic field B = Be z and an 
anisotropic harmonic potential. Choosing the gauge A = 
~B x r for the vector potential, the Hamiltonian can be 
expressed as H = H z + H xy , where the Hamiltonian 
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describes the z-direction with trap frequency u> z and 
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describes the x,y plane, where = */ w 2 + u 2 with oj; = 
^jfj. Equation @ can be derived from Eq. (0) without 
the 14 x tra term. Thus the eigenvalues of quantized levels 
for a boson are of the form 
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with n z = 0,1,2,..., n p = 0,1,2,..., and m = 0, ±1, 
±2,..., where some dimensionless variables are intro- 
duced: t n ^n p ,m =_e nz ,n p ,m/(fiwz), B = u>i/oj z and a = 
ljq/lo z . Note that £? is proportional to the external mag- 
netic field (synthetic magnetic field) for charged (neutral) 
particles. Generally there is no degeneration of energy 
levels, except that t he accide ntal dege neracy a rises when 
at least one of B, \J o? + B 2 and Bjy a 2 + B 2 is a ratio- 
nal number, eo = £o,o,o = 5 + V« 2 + B 2 is the ground 
state energy of one boson. 

Now we consider an assembly of N bosons, whose ther- 
modynamic potential is written as 51 = Six + ^0, where 
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with the dimensionless temperature T = ksT '/ '(Hu) z ) and 
chemical potential p — /j,/(hu} z ), and fio is the thermo- 
dynamic potentials for condensed bosons, fio is present 
only when there is a nonzero condensate described by a 
background field which reads 
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Here = KV^> + i^A^f? is the gauge-covariant deriva- 
tive and V(r) is the trapping potential as in Eqs. © 



and ([3]). The presence of a condensate is signalled by a 
nonzero value for ^ which satisfies 



i^- = — D 2 ^ - V(rW + n^ = 0. 
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<f can be expanded as *(r) = J2 n C nf n {r), where /„(r) 
are stationary state solutions to the Schrodinger equa- 
tion, 

- ^^ 2 /»(r) + V(r)/„(r) = E n f n (r). (9) 

and they form a complete set. For a ideal gas, the eigen- 
values E n coincide with the one particle energy described 
by Eq. ([5]). Combining Eqs. flS) and (HJ), the coefficients 
C n are determined by (E n — fi)C n — 0. The chemical po- 
tential is smaller than the ground state energy, /i < eo, 
above the BEC temperature T c , and therefore E n — u, 7^ 
and C n = for all levels. With the temperature decreas- 
ing, /x approaches eo until [i = eo at or below T c . Then the 
only solution C n= o 7^ exists at T < T c and the conden- 
sate wave function is nonzero correspondingly, *5f(r) = 
C fo(r), where /„ = 6j^(0, 1, a 2 p 2 )e-^ P 2 + a >- 2 ) is the 
normalized eigenfunction corresponding to the lowest en- 

ergy e . ^(a,j,x) = £ fc=0 n 7 ( 7 +i)...( 7 +fc-i) 1 ' 1S con " 
fluent hypergeometric function. 6 represents the normal- 
izing constant, a; = y/M\ui\/h and a z = \J Mco z /h. 

The BEC temperature and condensate fraction. - The 
particle number is derived from the thermodynamic po- 
tential by the standard procedure N = —dCl/dp, 
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Here E denotes the summation over all the eigen states 
but the ground state is not included in J2' ■ No is the 
number of particles condensed on the ground state, 
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Formally, the dimensionless chemical potential p is a 
function of the external magnetic field B and the tem- 
perature T, which can be obtained from Eq. (|TU)) . As 
discussed above, eo — p — > as T — > T+. Above T c , we 
obtain an analytical expression for the particle number 
by converting the summation in Eq. (I10[) into definite 
integral [25j . 
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where g 1 {z) is the polylogarithm function defined as 
g-y(z) = E^°=i e n i ■ ^ ne ser i es converge when 7 > 1 
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and z > 0. <7 7 (0) = £(7) is just the Riemann zeta func- 
tion. The BEC temperature T c is determined by setting 
e - p to be zero in Eq. JEJ), that is iV = T c 3 g 3 (0)/a 2 
where 53(0) = 1.202. So we have 
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This is just the BEC temperature for a Bose gas of TV par- 
ticles confined in an anisotropic trap with the frequencies 
cj in x,y plane and lj z in the z-direction. Since the one- 
particle energy spectra has already significantly changed 
by the external magnetic field, as suggested in Eq. ([5]), it 
is very surprising that the BEC temperature T c is irrele- 
vant to the magnetic field B. This result seems conflict 
with our intuition established in the study of supercon- 
ductivity, where the transition temperature decreases as 
the magnetic field is strengthened. Below T c , the number 
of thermal particles is Nt 
condensate fraction reads 



j^g 3 (0) = N(jr ) 3 , so the 
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Again, no varies as if the magnetic field does not exist. 

We now consider a Bose gas confined in a rotating 
frame, as described by Eqs. ((T|) and ©. Similarly, we 
obtain the BEC temperature, 
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With the rotation frequency ft rising up to u , the BEC 
temperature goes down to zero. This result is consistent 
with the predication that the ground state in the fast 
rotation regime undergoes a quantum transition to a non- 
condensed state as il — > ujo [3, EH . 

It is worth noting that Eq. (fl5|) is also helpful to un- 
derstand the magnetic properties of CBGs. With fi ap- 
proaching ujo , the trapping is canceled by the centrifugal 
effect and the rotating Bose gas reduces to a homogenous 
CBG in a constant external field, as Eq. (jTJ) indicates. 
It is well established that the BEC does no longer take 
place in this case pll - tzij l . 

The internal energy and specific heat. -The internal 
energy is calculated for two cases separately. For Case I 
that T > T c , the internal energy is given by 
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For Case II that T < T c , U is expressed as 
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FIG. 1: The specific heat C/N as a function of normalized 
temperature T/T c for charged Bose gases in magnetic field. 
This curve is universal to all the trapped Bose gases. 



Correspondingly, the specific heat can be derived from 
the internal energy, 
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(18) 

At T > T c , e — fJ. satisfies Eq. (fT2|) and it is just the 
function of temperature T, while eo — p, — at T < T c . 
In both cases, eo — p is independent on the magnetic field 
B. So, the S-dependency of U in Eqs. {M} and (JTJJ) 
is just due to the eo-^ term where eo is in relation to B. 
Nevertheless, since the eo-/V term remains constant as the 
temperature varies, the specific heat has no relation with 
B. Figure [TJ shows the specific heat as a function of the 
normalized temperature, which displays a discontinuity 
at the BEC temperature. 

The Landau diamagnetism. - We now calculate the 
magnetization according to M = —dVl/dB, which yields 
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at T > T c . Below T c , M consists of two parts. One is 
due to the condensed particles, 
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The other is the contribution from thermal particles 
dn T N T B 



(20) 



M T = 



dB Va 2 + B 2 
Note that the summation of both parts, 

NB 

M = Mn + M T 
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(22) 
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just amounts to the result above T c . In addition, the 
susceptibility has the following form, 
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(23) 



It needs points out that both the magnetization and sus- 
ceptibility keep invariant at all temperatures. According 
to Eqs. (|19p . (|2"2")l and (|2"3"|) . the magnetization and sus- 
ceptibility are negative, reflecting that the gas exhibits 
the Landau diamagnetism. 

The diamagnetization for a homogeneous CBG in the 
ma gne tic field has been investigated by several groups 
[22H24| . The magnetization M varies with the tempera- 
ture and the diamagnetism is stronger at lower temper- 
atures. Especially, below the BEC temperature, M does 
not vanish as the external field B is reduced to zero, 
which is the evidence of the Meissner-Ochsenfeld effect. 
Nevertheless, present results shows that M vanishes as 
B — > at all temperatures, implying that the Meissner- 
Ochsenfeld effect might not exist in a trapped CBG. The 
trapping potential brings about significant changes to the 
physics of the charged Bose gases. 

In summary, we show that the charged Bose gas con- 
fined in a harmonic trap undergoes Bose-Einstein con- 
densation at a critical temperature determined by the 
trapping potentials and irrelevant to the magnetic field, 
although the application of the field changes the energy 



spectrum of bosons. The specific heat is also indepen- 
dent of the external field. Moreover, we find that the 
Landau diamagnetization is strengthened with the mag- 
netic field, but keeps unchanging at all temperatures. 
Our results are applicable to the neutral atomic Bose gas 
which is caused to rotate either by a rotating frame or due 
to a synthetic magnetic field. The latter behaves like a 
charged Bose gas in magnetic field, while the Bose gas in 
a rotating frame displays distinct thermodynamic behav- 
iors. Its Bose-Einstein condensation temperature does 
no longer hold a constant but decreases with the rotat- 
ing frequency increasing. The condensation does not oc- 
cur once the rotation frequency approaches the trap fre- 
quency, consistent with the predication that the ground 
state in the fast rotation regime is a non-condensed state. 

At last, we briefly discuss the possibility to realize 
charged Bose gases in experiments. It is already possible 
to create ultracold plasmas by photoionization of laser- 
cooled neutron atoms [26| . The temperatures of electrons 
and ions are as low as 100 mK and 10 /iK, respectively. 
The ions can be regarded as charged bosons if their spin 
is an integer. 
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